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Instabilities of Non-Abelian Vortices in Dense QCD 
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We construct a low-energy effective theory describing non-Abelian vortices in the color super- 
conducting quark matter under stress. We demonstrate that all the vortices are radically unstable 
against decay into the only one type of vortices due to the potential term induced by the explicit 
flavor symmetry breaking by the strange quark mass. A simple analytical estimate for the lifetime 
of unstable vortices is provided under the controlled weak-coupling calculations. We briefly discuss 
the (non)existence of magnetic monopoles at high density. 
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Introduction. — Topological or quantized vortices com- 
monly arise in a wide area of physics from condensed 
matter physics and cosmology to particle physics In 
the context of nuclear physics, the dynamical breaking 
of U(1)b baryon number due to the neutron superflu- 
idity in nuclear matter gives rise to topologically sta- 
ble vortices characterized by the first homotopy group 
7ri[U(l)s] — Z. They are phenomenologically impor- 
tant since the sudden increase of the rotation of neutron 
stars, the so-called glitches, may be attributed to the un- 
pinning of vortices which releases an angular momentum 
transfer from the nuclear "mantle" to the outer crust Q . 
Topological superfluid vortices also emerge in the color 
superconducting quark matter Q presumably existing in 
the "core" of neutron stars: the U(1)b symmetry is bro- 
ken by the condensation of diquark pairs in the color- 
flavor locked (CFL) phase ^] which is indeed shown to be 
the most stable ground state at asymptotic high density 
in quantum chromodynamics (QCD). Recently, however, 
it has been found that minimal topological vortices in 
quark matter are not U(1)b vortices [5] but non-Abelian 
vortices @ referred to as the semisuperfluid vortices [?[, 
which have only winding number 1/3 inside U(1)b. Ac- 
tually it is energetically favorable for a single U(1)b vor- 
tex to split into three (red, green, and blue) non-Abelian 
vortices. At first glance, one may expect that all the 
resultant three non-Abelian vortices are stable. 

In this Letter, we show that these remaining non- 
Abelian vortices are still unstable against decay into the 
only one type of stable vortices when the effect of nonzero 
strange quark mass rUg is taken into account. In order 
to elucidate the (in)stabilities of non-Abelian vortices 
in a model-independent manner, we use the Ginzburg- 
Landau (GL) approach near the transition temperature 
Tc, and construct the low-energy effective theory of vor- 
tices with the potential term induced by the explicit 
breaking of flavor symmetry. Owing to the asymptotic 
freedom of QCD, all the calculations throughout this Let- 
ter are under theoretical control at high density regime 



where the QCD coupling constant is weak. We remark 
that the existence of non-Abelian vortices by itself does 
not rely on the domain of applicability of the GL La- 
grangian, but only on the dynamical symmetry breaking 
induced by the diquark condensation. On the other hand, 
the symmetry argument is not enough to ensure their sta- 
bilities which depend on the details of the dynamics. In 
the following, we neglect the effect of U(1)em electromag- 
netism since the mixing between broken SU(3)c color 
and U(1)em is sufficiently small at high density. The 
generalization to include the effect is straightforward. 

Ginzburg- Landau Lagrangian. — We consider the di- 
quark pairing in the most attractive CFL and spin- 
parity 0+ channel [i]: {^hYa ^ eabc^ijk{(qL)iC{qL)'^) and 

('J'fl)a ^ <^ahc^tjk{{qR)iC{qR)'^^), where i,j,k {a,b,c) are 
flavor (color) indices and C is the charge conjugation op- 
erator. Here we take — —^r = $ so that the ground 
state is the positive parity state. 

The time-dependent Ginzburg-Landau (TDGL) La- 
grangian up to the second order in time and space deriva- 
tives [8] at large quark chemical potential /i 3> rus ^ 
rriu.d — near Tc is given by [jl. [lo|: 



-GL 



+Tr {Kb^^Vo^ + H.c.) - ^F^.F^"" ~ Fgl, 



Tr 



+/3i [Tr($t$)]%/32Tr [($^$)2] , 



(1) 



where X'^j^ 



igsAy,^, F^^ 
= idiag(0,l,- 



- d^Au — dijA^ — 
and = ^diag(0, 1, — 1). JCd is a dis- 
sipative term reflecting the decay of Cooper pairs into 
fermionic excitations. The e terms originate from a Fermi 
surface splitting due to the nonzero strange quark mass 
together with the constraints of electric and color charge 
neutrality and weak interaction equilibration [lo| . 

The GL parameters Kq,^, a, /3i,2, and e are obtained 



2 



from the weak-coupling calculations [13] : 

3 12(7rTc)2 ^^V^ ^Tc'^^ 

where N{fj,) = fj? /{2t:'^) is the density of state at the 
Fermi surface andTc = 2^/'^ei l^/ TT is the critical temper- 
ature of the CFL phase in the absence of mg. Kq and 
Kd have not been calculated in the literature, but can 



be derived following the same procedure of Ref. [U 



The ground state of the GL potential is given by $ = 

\ -|l/2 

~^ ~ W j - ^^sj = diag(Ai, A2, A3), where 
the gap parameters Ai, A2, and A3 denote down- 
strange, strange-up, and up-down Cooper pairs, respec- 
tively. Due to the gap orderi ng, A3 > Ai > A2, the 
symmetry breaking pattern is [lOj 



SU(3)c X SU(3)L,fl X U(l)s ^ SU(3)c^ 



U(l)|..(3) 



For clarity and completeness, we will first neglect the £^3 
term and later treat it as a perturbation. Without the 
6X3 term, the order parameter is given by $ = AI3 = 



8/3 12^ 

Mass spectra around this ground state are 



m2j = 2.g2A2if3 



2e 

T 



4^A2 



(4) 



where tuq is the mass of the gluons, mi and ms are 
the masses of quarks in the 3 (g) 3 = 1 © 8 representa- 
tion under the unbroken SU(3)c+_f symmetry, respec- 
tively. From Eqs. ([2]) and ([4]), we have tog ~ ffsA* and 
TOi ~ 2to8 ~ A; then, the relation g^/i 3> A at high 
density indicates that the CFL phase is a typc-I super- 
conductor 1^. Note that non-Abelian vortices can ap- 



pear even in this type-I system, since their interactions 
are repulsive at large distances due to the exchange of 
the Nambu-Goldstone (NG) boson associated with the 
U(l)s symmetry breaking flS^. This is in contrast to 
the case of the metallic (Abelian) superconductor where 
vortices can exist only in a type-II system. Non-Abelian 
vortices are rather superfluid vortices; they are created 
under a rapid rotation. 

Profiles of non-Abelian vortices. — Corresponding to 
the three types of vortices, the order parameter $ asymp- 
totically behaves as 



diag(0,*,*) 
$ ' ^"> { diag(*,0, *) 
diag(*,*,0) 



diag(e*^Ai,A2,A3) 
$1^<| diag Ai,e*^A2,A3 
diag(Ai,A2,e^^A3) 



where (r, 6) is the polar coordinate and "*" stand for 
some nonzero constants. All the above three asymp- 
totic forms at infinity can be brought into a unique 



form $ e*^/3 diag(Ai, A2, A3) by regular SU(3)c 

gauge transformations |13|. The overall phase e*^/^ man- 
ifestly shows that the non-Abelian vortex winds 27r/3 
inside U(1)b. Their tensions logarithmically diverge as 
T ~ log ^ -h C(l) where L is a long-distance cutoff 
and ro is a short-distance cutoff. 

Let us take a diagonal ansatz for a single vortex 



A e 



Fjr) 
V3 



To 



(5) 
(6) 

-1), and 
(3,0,0) 



with To = ^diag(l,l,l), T3 = -i=diag(0,l, 
Tg = ^diag(-2,l,l). We impose {F,G,h) - 
as r — 00 to satisfy $ — >■ AI3. The single-valuedness 
condition for $ requires {v^i^vs) = (0, 1), ~\)- 

In the presence of each vortex, the remaining 
SU(3)c+_F symmetry is further broken down to [U(l) x 
S\]{2)]c+F ■ Hence, the vortex solution is labeled by the 
NG modes (or the orientational modes) living on the 
coset space SU(3)c+f/[U(1) x S\](2)]c+f ^ CP^, which 



we parametrize by introducing 
1, (j) ~ e*"0) defined as U 



C/t = 



^13. The most general solution can be obtained by 



acting U G SU(3)c+_f as <i> 
"F(r 



Ae^ 



V3 



1 f T-J 
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U^W and A. 
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UAiW: 
(7) 
(8) 



(0,1) 



For concreteness, let us choose (1/3, i/g) 
reference solution. Equations of motion for the profile 
functions read 



9 +^ 
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9_ 
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1t!_ 

r 



TO? 



[g^{h-l) + f{2h+l)\^Q, 



(9) 



with / = \{F + 2G) and 5 = \{F - G). These equa- 
tions are solved with the boundary conditions, (/, g, h) — > 
(1, 1, 0) as r -J> 00 and (/, g', h) (0, 0, 1) as r ^ 0. 

Low-energy effective theory. — The NG modes 4> G CP^ 
propagate along the non-Abelian vortex string. The 
philosophy of constructing the low-energy effective La- 
grangian is similar to that of the chiral perturbation 
theory (ChPT) describing the low-energy dynamics of 
QCD. Remembering that the ChPT is constrained by 
the [SU(7V/)l X SU(iV/)fl]/SU(7V/) symmetry, the form 
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of the Lagrangian in our case is determined solely by the 
SU(3)/[U(1) X SU(2)] symmetry, and is described by the 
CP^ nonlinear sigma model |15| : 

a=0,3 

where is promoted to a dynamical field as (/) — ?> 4>{xo, xs) 
depending on the vortex world-sheet coordinates xq and 
X3. Ka are the stiffness parameters in Eq. ([1]) and we 
have only one unknown constant C. Note that the K^:,- 
term in Eq. ([1]) gives no contribution to Eq. ()10l) since it 
is traceless in the vortex background solutions [15]. 

In order to determine the constant C, we have to go 
back to the original GL Lagrangian ((TJ and we have to 
know the ip dependences of $ and . It is easy for $ and 
■^1=1,2 because we have already solved background vor- 
tex solutions as ^{xi^2', (pixo.s)) and Ai^2ixi^2T 'Pi^o.a))- 
The missing part is Ao,3(0(a;o,3)) which vanishes in the 
background solutions. Therefore we make an ansatz in an 
appropriate gauge following Ref. ^Q(a;i, 2; 0(2:0,3)) = 
^^^[00^, 9q,(00^)] (Qf = 0, 3) where p{r) is an unknown 
function. Then we finally arrive at 



C = 



An [ r 
^7 '^2 



4((l-p)(/-.9)' + y(/' + 5^) 



(11) 



0, one finds 



where p should be determined so that the integral pT|) 
is minimized. Using the Euler-Lagrange equation for p, 

+ (l-p)^-^ [U^+g^)p-{f-gf] 

that C is indeed finite and (j) is normalizable |14| . 

Unstable non-Abelian vortices. — We now turn on the 
eX^ term and consider the regime e ^ a, which allows 
for an analytical treatment. Since this term explicitly 
breaks SU(3)c+f symmetry, the NG modes in Eq. 
are lifted via an effective potential over the CP^ space. 
Let us consider a single vortex whose field configuration 
satisfies Eq. ([9]) . Variations of its tension can be thought 
of as the potential 

^ejcPx Tr [$^^3$] = i^(|03|' - 1021'), (12) 

where we have used |0iP + |02p + |03p = 1 and have 
defined 



D = TreA^ / dr r{g^ - f). 



(13) 



Note that D is positive and finite because g — f is al- 
ways positive and gets exponentially small as going away 
from the vortex jlj|; thus, the effective potential is well- 
defined. 

The effective potential in CP^ space is shown in Fig.[T] 
Since the potential has one minimum at (0i, 02,03) = 



1.0 

0.8 
. 0.6 
" 0.4 

0.2 



0.0 




I 



0.0 0.2 0.4 0.6 0.8 1.0 



+1.0D 

+0.8 D 
+0.6 D 
+0.4 D 
+0.2 D 


-0.2 D 
-0.4 D 
-0.6 D 
-0.8 D 
-l.OD 



FIG. 1: (Color) Contour plot of the effective potential for the 



CP-" NG modes in the |02|' 
the height of the potential. 



plane. The color represents 



(0,1,0), any vortices away from (0,1,0) are unstable 
against decay into the (0, 1,0) vortex. This matches the 
fact that the pairing gap A2 is smaller than Ai and A3 
so that the vortex whose string tension is proportional 
to A2 is easier to be created than others; the details of 
the dynamics even suggest that the (1, 0, 0) and (0, 0, 1) 
vortices are no longer local minima. 

Let us estimate the lifetime of unstable vortices. As an 
example, we consider the decay from the (1,0,0) vortex 
at the left-bottom corner of Fig. [T]to the (0, 1, 0) vortex 
at the right-bottom corner. The discussion here holds 
for the (0, 0, 1) vortex. In what follows, we set 03 = 0, 
implying that we will consider a CP^ submanifold (cor- 
responding to the bottom edge of Fig. [1]) inside CP^. 
It is useful to introduce an inhomogeneous coordinate 
u{t) e CPi by (01,02) = il/yTTW,u/y^T+W)- 
Then the low-energy effective Lagrangian can be rewrit- 
ten as 



•^cpi — CKq 



(1 + H2)2 



D- 



1 



(14) 



A typical time scale of this equation is r = ^/CKq/D. 

In principle, we can numerically calculate r for each 
p. Here we provide a simple analytical estimate instead. 
Since the profile function / {g, h, and p) increases (de- 
crease) with a typical scale r ^ A^^ for tog ^ Tni,s 
14j . we find C ^ {p/^Y from Eq. (fTTj) . Furthermore, 
D is estimated from Eq. p3)) as D ^ e ~ m^log(/i/A). 
Therefore the lifetime of unstable vortices is given by 



TOs ^ri{p/A), ri{x) = x"^ (log xY 



(15) 



In the limit — ^ 0, r — >■ cx) as anticipated. 

(Non) existence of magnetic monopoles. — Let us dis- 
cuss the (non)existence of magnetic monopoles in QCD at 
high density. One may expect that the symmetry break- 
ing pattern ([3]) would support the magnetic monopoles 
characterized by 7r2[SU(3)/U(l)^] —1?. If so, monopoles 
must be confined due to the color Meissner effect of 
the color superconductivity because it is in the Higgs 
phase. In fact, such a confined monopole exists in the 
N = 2 supersymmetric QCD in the Higgs phase with 
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the same symmetry breaking pattern ^ , where mag- 
netic fluxes are squeezed into vortex strings confining the 
monopole from both sides. This composite object has 
been understood as a kink in the low-energy effective 
world-sheet theory on the vortex string with a suitable 
potential term admitting more than or equal to two min- 
ima. If the low-energy theory (ITOl) in our case had a 
potential similar to supersymmetric QCD, this would re- 
alize the dual of the confinement scenario advocated in 
the QCD vacuum where monopoles are condensed and 
quarks are confined 18[. However, this is not the case. 
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The potential (112)) has only one minimum and allows no 
kink solutions but implies the instabilities of non-Abelian 
vortices instead, as we have seen. 

One should note that this conclusion may not be 
valid if one includes the nonperturbative quantum effects 
which account for the mass gap of NG modes as indi- 
cated by the Coleman-Mermin- Wagner theorem in two 
dimensions. Actually, such effects may lead to multiple 
local minima in the potential, and thus, the monopole- 
antimonopole meson attached to the vortex 1^ 13 • 
the original four-dimensional GL theory at sufficiently 
high density, instanton effects are highly suppressed due 
to the screening of instantons together with the asymp- 
totic freedom of QCD [io^, and another mechanism re- 
sponsible for the quantum effects should be present. We 
will defer this issue to a future work. 

Discussion. — It is interesting to investigate possible as- 
trophysical implications of our results. When the core 
of a neutron star cools down below the critical tempera- 
ture of the CFL phase, a network of non-Abelian vortices 
will be formed by the Kibble mechanism. Remarkably, 
the extrapolation of our formula (jl5p to the intermedi- 
ate density regime relevant to the core of neutron stars 
in ~ 500 MeV) with A - 10 MeV and rus ~ 150 MeV 
suggests that all the vortices decay radically with the 
lifetime of order r ~ 10~^^ second. Although this result 
should be taken with some care due to the uncertainty 
of numerical factor in Eq. (|15p , it is reasonable to expect 
that only one type of non-Abelian vortices, which cor- 
respond to the point (0, 1,0) in the CP^ space, survive 
as a response to the rotation of neutron stars in reality. 
The other decaying non-Abelian vortices will emit NG 
bosons, quarks, gluons, or photons during thermal evo- 
lution of neutron stars. In relation to the glitch phenom- 
ena, it would be also important to understand how the 
Abelian U(l)s vortices in hadronic matter are connected 
to the stable non-Abelian vortices in color superconduct- 
ing quark matter in the interior of neutron stars. This 
may be relevant to the question of continuity of hadronic 
matter and quark matter 2l|, |22| . 
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